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ABSTRACT 

We studied the X-ray timing and spectral variability of the X-ray source Sw J1644+57, 
a candidate for a tidal disruption event. We have separated the long-term trend (an 
initial decline followed by a plateau) from the short-term dips in the Swift light-curve. 
Power spectra and Lomb-Scarglc pcriodograms hint at possible periodic modulation. 
By using structure function analysis, we have shown that the dips were not random but 
occurred preferentially at time intervals ~ (2.3,4.5,9) x 10 5 s and their higher-order 
multiples. After the plateau epoch, dipping resumed at « (0.7, 1.4) x 10 6 s and their 
multiples. We have also found that the X-ray spectrum became much softer during 
each of the early dip, while the spectrum outside the dips became mildly harder in its 
long-term evolution. We propose that the jet in the system undergoes precession and 
nutation, which causes the collimated core of the jet briefly to go out of our line of 
sight. The combined effects of precession and nutation provide a natural explanation 
for the peculiar patterns of the dips. We interpret the slow hardening of the baseline 
flux as a transition from an extended, optically thin emission region to a compact, 
more opaque emission core at the base of the jet. 

Key words: accretion, accretion discs — galaxies: jets — methods: data analysis — 
galaxies: active — X-rays: individual: Sw J1644+57 — black hole physics. 



1 INTRODUCTION 

Swift J164449. 3+573451 (henceforth, Sw J1644+57) was 
discovered by the Swift Burst Alert T elescope (BAT; 
iGehrels et al.ll2004 [Barthelmv et aJ1l2005l) on 2 011 March 
28 (|Cummings et al.ll201ll ; IBurrows et al.l 120111 ). The long 
duration of the X-ray emission and flaring events (still on- 
going after 8 months), and the spatial coincidence with the 
nucleus of a galaxy at redshift z = 0.354 (l uminosity dis- 
tance « 5.7 x 10 27 cm: iLevan et aUl2011bl lal; iThoene et all 
l201ll ) made it clear that it was not a Gamma-Ray-Burst, 
and was instead associated to some kind of sud den accretion 
event onto a supermassive black hole (BH) jBloom et alj 
l2011bl ; IBurrows et "ail 1201 ll ; ILevan et al.ll2011ah . In partic- 
ular, the precise position from Very Large Array and Very 
Long Baseline Array radio observations clinched the identi- 
fication of the X-ray source with the nucleus of the host 
galaxy (|Zauderer et al.l l201ll ). The lack of previous evi- 
dence of nuclear activity (AGN) from the same source, 
and the short timescale fo r the outburst rise (a few days: 
iKrimm fc Barthelmv! 1201 ll ) are difficult to reconcile with 
changes in the large-scale accretion flow through an ac- 
cretion disc, but this scenario is not completely ruled 
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out (jBurrows et alj l201ll ) . The most favoured scenario is 
a tidal disrupti on event deep insid e the sphere of influ- 



ence of the BH JBloom et al 



Ca 



Troia fc Lodatol l20 



201 lallbl; IBurrows et al.ll201ll : 



Socratesll201ll ). The nature 



of the star being disrupted (and therefore the radius at which 
the event occurs and the timescale on which the stellar ma- 
terial circularised into a disc and accreted) are still matters 
of intense debate. For example, it was argued that if the dis- 
rupted star were a white dwarf, the smaller tidal disruption 
radius and characteristic timescale would be more consis- 
tent w ith the observed duration of the flares (|Krolik fc Piranl 

Eon]). 

There are no direct measurements of the BH mass 
(M.). Indirect order-of-magnitude constraints consist of: 
an upper limit to the BH mass that permits tidal dis- 



ruption events outside the event horizons (|Rees! 1 19881 ); 
a minimum variability timescale ~ 100 s that sets an 
3er limit to the li g ht-cro ssing timescale (|Bloom et alj 



uppe 



l2011bl ; IBurrows et al.l l201ll ); and well known empirical 
relations between M, and the host galaxian environ- 
ment (whatever causes those trends, e.g. ISilk fc Reeslll998l ; 
iJahnke fc Macciol l201ll) . For exampl e, based on t he op - 
tical luminosity of the host galaxy, ILevan et al.l (|2011al ) 
estimate a spheroidal stellar mass M sph ~ 10 9 -10 10 M Q . 
Us ing the log-linear sp heroidal mass-BH mass relation 
of iBennert et all (|201lf ) (see also iMagorrian et al.l Il998l ; 
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lLauer et alj|2007l : iKormendv. Bender" fc CornellllibTlh , this 
implies a likely BH mass 2 x 10 6 M Q < M BH < 10 7 M Q 
ijLevan et a,l.ll2011ah ; adop ting instead the bent M sp h-M] 



relation of 
BH masses 



Graham! (|2012h . it allows also a lower ranee of 



10 s M Q < Mbh < 10 7 M Q . The latter range is 
more consiste nt with the BH mass log (M. /Mp>) — 5.5 ± 1.1 
estimated bv lMiller fe GiiltekinI j 201 ll ). based on empirical 
"fundamental plane" relations between radio and X-ray lu- 
minosities of accreting BHs. 

The relatively small BH mass implies a bolometric Ed- 
dington luminosity < 10 45 erg s -1 , and a 0.3-10 keV lu- 
minosity a few times smaller. The average isotropic lumi- 
nosity in the 0.3-10 keV band was « 10 47 erg s" 1 in the 
first few weeks after the initial outburst; several months 
(~ 10 7 s) later, it was still ~ a few xlO 45 erg s -1 , well 
above the Eddington limit. This is s trong evidence that 
the emission is beamed towards us |Bloom et all l2011bl ; 
iBurrows et"aH l201ll ; |I evan et al.ll2011ai r and is probably 
associated with a relativistic jet, like in blazars. A rela- 
tivistic jet has also been invoked to explain the radio tran- 
sient event, interpreted as an external shock in the gas 
surrounding Sw J1644+57 (uBower Bloom fc Cenkol 1201 ll ; 
IZauderer etai1l201ll ; lBurrows et aljhoill ). The Lorentz fac- 
tor of the jet is highly uncertain , with estimates varying from 
r » 2 - 5 IZauderer et aTll201lf) to F w 10-20 (|Bloom et all 
l2011bl ; IBurrows et al.l l201lh . The rarity of the event, with 
one such outburslQ discovered in 7 years of Swift operations, 
is consistent with the theoretical expectations of tidal dis- 
ruption rates and jet opening angle 6 ~ 5°, corresponding 
to T « 1/0 f« 10. 

The origin of the photon emission is still unclear. There 
are clearly at least two prominent peaks in the spectral en- 
ergy distri bution: one in the fa r IR and one in the h ard X- 
ray band (|Bloom et al.ll2011bl ; IBurrows et all^Onl ). They 
can be modelled as direct synchrotron emission (single- 
component model) from radio to X-rays, with strong dust 
extinction in the optical/UV band. Alternatively, the ra- 
dio/IR peak is the direct synchrotron emission and the X- 
ray peak is due to inverse Compton scattering of external 
photons, most likely disc photons (two-component blazar 
model). A third possibility is that the X-ray emission is due 
to inverse Compton emission at the base of the jet, while 
the radio/IR synchrotron emission comes from the forward 
shock at the interfac e between the head of the jet and the 
interstellar medium (|Bloom et al.ll2011bl ). 

The X-ray spectrum is well fitted by a simple absorbed 
power-l aw, although more complex models were also dis- 
cussed (|Burrows et alj201ll ) . The avera ge photon index dur- 



ing the first » 2 x 10 s s is T « 1.8 (jLevan et al.l l2011al ; 



IBurrows et al.ll201ll ). However, the physical meaning of this 
average value has to be interpreted more carefully, as there 
are strong variations in hardness and photon in dex between 
flares, with F changing between m 1.3 and 3 (|Levan et al.l 
l2011al ). In particul ar, the photon index is harder when the 



source is brighter l|Kennea et al.ll201ll; IBloom et aT1l2011al ; 
iLevan et afll2011al ; uBurrows et al.ll201lh . 

In the Swift X-Ray Telescope (XRT; IBurrows et alj 



10~ 7 -10" 4 Hz frequency 



120051 ) band, 0.3-10 keV, the declin e in flux was ini- 
tially consistent wit h the t~ 5//3 decay |Levan et al.l l2011al : 
IBloom et al.l l2011bT ) expected for the spreading and ac- 
cretion of fa l lback m aterial after a tidal disruption event. 
Bloom et all j201 lbl ) used the early Swift/XRT observations 
(up to t ~ 4.5 x 10 6 s) to study the power spectral density 
over the 0.1-100 mHz frequency range: they found no signif- 
icant feature that is not associated with the orbital period of 
the spacecraft. By contrast, a series of XMM-Newton Target- 
of-Opportunity observations during the initial decline phase 
suggested the pre sence of a 4.7-mHz (aa 200 s) quasi pe- 
riodic oscillation (|Miller fc Strohmaverl 120111 ) . and showed 
short-term variability over the 0.4-100 mHz range. If the 
putative quasi-periodic oscillation corresponds to the Kep- 
lerian frequency at the innermost stable circular orbit of a 
Schwarzschild BH, it would imply a BH mass ~5x 1O 5 M0 
|Miller fc Strohmaverlkoill '). 

Possible variability on the 
range has not been investigated in much detail yet. The on- 
going series of Swift/XRT monitoring observations, which 
have now followed the source for « 2 x 10 7 s, offer 
the best chance for this study. Luckily, eight months af- 
ter the outburst, the source is still detected at a flux ~ 
10~ n erg cm~ 2 s _1 in the 0.3-10 keV band, declining only 
very slowly. In this paper, we examine the X-ray variability 
over the 10 4 -10 7 s timescale, looking for phenomenological 
patterns (characteristic variability timescales, spectral evo- 
lution) that can help us test physical scenarios. In particular, 
we want to determine whether the dips that characterise the 
X-ray light-curve are random or have some periodicity, and 
we discuss the X-ray changes during the dips. 



2 DATA ANALYSIS 

Sw J1644+57 has been monitored by the XRT several times 
a day, every day since 2011 March 28. We used the on- 
line S wift/XRT data product generatoi0 |Evans et al.ll2007i . 
120091 ) to extract light-curves in different energy bands, and 
spectra (including background and ancillary response files) 
in different time intervals; we selected grade 0-12 events. We 
downloaded the suitable spectral response files for single and 
double events in photon-counting (PC) mode and window- 
timing (WT) mode from t he latest Cali bration Database. 
We used XSPEC Version 12 ( Arnaudlll996l') for spectral anal- 
ysis, and standard FTOOLiJ 3 |Blackburnlll995h tasks for pre- 
liminary timing analysis — for example, for defining the time 
intervals that were used to extract intensity-selected spectra. 
For more advanced timing analyses, we used a Lomb-Scargle 
periodogram and structure function analysis, (discussed in 
Section 3). 

A caveat for timing analysis is that the observational 
data are irregularly sampled; the duration of each snap- 
shot observation in either WT or PC mode (typically, a few 
hundred seconds) and the temporal gaps between observa- 
tions are not constant. Irregular, gappy data sequences are 
accommodated within the formulation of LS-periodograms 
(£| |3.2.2[) . where data are at discrete events rather than long 



1 It has been argued that there may be a sec ond object discov- 
ered by Swift in this class l lCenko e t al.ll201ll). and another one 
discovered by XMM-Newton iLin et al.ll201ll 'l. 



2 Including the new treatment of the vignetting correction, in- 
troduced after 2011 August 5. 

3 http: / /heasarc . gsf c .nasa. gov/f tools 
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Figure 1. Top panel: Swift/XKX 0.3-10 kcV light-curve between 
2011 March 28 and 2011 November 28, plotted on a linear time 
scale. Datapoints in red are taken in PC mode; those in magenta 
are in WT mode. The dashed curve is a phcnomcnological func- 
tion used to de-trend the data and "normalise" the continuum. 
Bottom panel: same as in the top panel, but on a logarithmic 
time scale. 



temporal bins. Such imperfections however are unnatural to 
conventional Fourier (§ 13.2. ip and structure function meth- 
ods (§ 13.2.31) . There are various strategies that can circum- 
vent this issue: we considered the hyphen method, the zigzag 
method and the trapezoid method. In the hyphen method 
the gaps are practically omitted. The flux levels are set to 
be locally constant for the durations of the data bins. In the 
zigzag method, a data point is set at the mid-time of each 
observational bin, and we connect each point linearly to the 
next consecutive data point. The numerical light curve is 
piecewise linear and resembles a polygonal landscape. In a 
variant trapezoid method, the observational bins are treated 
as piecewise flat segments as in the hyphen method but the 
gaps are connected directly with diagonal lines. 

As the covering factor of the bins is small in the hyphen 
method, unless the data are numerically bridged, we found 
that the calculated power spectra and structure functions 
tend to blur any useful information in the noise. Instead, 
we found that the zigzag method tends to provide robust 
results for power spectra and structure functions. The results 




Frequency (Hz) 

Figure 2. Power spectrum of the Swift/XRT light curve (com- 
bining data in WT and PC mode), excluding the first 2 X 10 5 s. 
The input data have be en binned to 5000 s intervals. The power 
spectrum is normalised jLeahv et alj|l983h such that its integral 
gives the squared rms fractional variability (i.e., the Y axis is in 
units of (rms) 2 Hz -1 ), and the expected white noise level Rj 2. 
Two features at /iHz frequencies (labelled in the plot) arc sugges- 
tive of periodicities on timescales fa (2.3 ±0.3) X 10 s s and about 
twice that value. 

obtained with the trapezoid method are indistinguishable 
from those of the zigzag method. Thus, we present here only 
the results from analyses done with the zigzag method. 



3 MAIN RESULTS 

3.1 Different phases in the lightcurve 

Figure [T] shows the XRT 0.3-10 keV light-curve, binned by 
snapshot (typical exposure duration ~ a few 100 s). The 
luminosity evolution shows a series of phases with different 
phenomenological properties. During the first three months 
after the outburst — in particular, at times 2 x 10 5 s < t < 
9 x 10 6 s — the baseline trend is adequately fitted by an ex- 
ponential decay, with an e-folding timescale of « 2.0 x 10 6 s. 
Flares and co mplex dips are sup erposed onto this trend. 
In other work l|Levan et al.ll2011al ). a more canonical tidal- 
disruption decline scaling ~ t~ 5 ^ 3 was fitted to the same 
section of the light-curve; the difference is mostly due to 
a different definition of the baseline. Later we will show a 
change in the variability at t ~ 4.7 x 10 6 s, which distin- 
guishes "early" and "late" epochs of the declining stage. For 
t > 9 x 10 6 s, the decay stopped and the source appeared to 
settle on a "plateau" (XRT count rate « 0.1 ct s _1 ), with- 
out any dips. Broader dips resumed for t > 12.5 x 10 6 s (the 
"recent" epoch), but the baseline level has not significantly 
declined below ~ 0.1 ct s _1 . 

Because of the different behaviour of the light-curve at 
different epochs, we shall investigate the short-term variabil- 
ity separately in the different phases. We also note that the 
state transitions are each a continuous, gradual evolution 
of flux, timing and spectral properties over time. Therefore 
the choice of start and end times defining the different epoch 
sub-intervals are somewhat arbitrary (to within a few 10 s) 
and do not affect the main conclusions of our study. 
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Figure 3. Lomb-Scarglc pcriodograms of the light-curve in four time intervals of interest. The data have not been de-trended. In the 
first, second and fourth of these intervals, there are plausible periodicities: significant at the 1% FAP level (with two definitions of this 
threshold marked as horizontal lines). 



3.2 Search for periodicities 

3.2.1 Fourier analysis 

Fourier techniques, as implemented for example in the 
XRO NOS powspec timing analysis task (|Stella fe Angelina 
are widely used for the timing analysis of X-ray 
lightcurves in AGN and X-ray binaries. The standard tech- 
nique of dividing the lightcurve into multiple intervals, cal- 
culating the power density spectrum in each interval, and 
taking the mean of the power density values (and cor- 
responding standa rd deviations) i s applicable on ly if the 
process is erg odic <|Priestlev|[l"988l ; lGuidorzill201ll ), so that 
time averages can be substituted for ensemble averages. In 
the case of highly non-stationary events like Sw J1644+57, 
Gamma-ray bursts, or X-ray flares, it does not make sense to 
subdivide the light-curve into sub-intervals, and the power 
density spectrum has to be calculated o ver the entire dura- 
tion of the observation l|Guidorzill201ll ). We did that, and 
then calculated the standard deviation for the power den- 
sity at each fre quency bin with the procedure outlined in 
iGuidorzil (|201ll ). 

Figure [2] shows the power spectrum of the 0.3-10 keV 
light-curve, excluding the first 2 x 10 5 s (that is, after the 
huge initial flares have subsided). The root-mean-square 



(rms) power rises at lower frequencies, as ~ v^ 1 . This is 
mostly due to the long-term dimming trend. There appear 
to be two unresolved features at frequencies ~ 2.2/iHz and 
4.5/iHz; however, their power is comparable to the statis- 
tical uncertainty. This prevents us from firmly concluding 
whether there are true underlying periodicities or whether 
those features were just statistical noise. Moreover, Fourier 
analysis techniques are better suited to time series that are 
cquispaced and with no time gaps. Neither condition is true 
in our case. The complex shape of the window and sampling 
functions may introduce side lobes (spectral leakage) in the 
discrete Fourier transform. In summary, we consider the fea- 
tures observed in the power density spectrum as suggestive 
of possible periodic signals, but not yet statistically signifi- 
cant, because of the shortcomings of the Fourier technique. 



3.2.2 Lomb-Scargle periodogram 

As an imp r ovem ent o ver Four i er pow er spectrum analy- 
ses, iLombl (1 19761 ) and IScargld (|l982l ) introduced a peri- 
odogram technique that proves to be more robust in the 
detection of periodicitie s in irregularly sampled light curves. 
iPress fc Rvbickil (|l989T l accelerated the method, with a nu- 
merical modification based on fast Fourier transforms. The 
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Lomb-Scargle periodogram is equivalent to a best-fit anal- 
ysis of the data with a single sinusoidal function. The 
phase and amplitude emerge directly for the implied fit. 
The technique is observationally applied to diverse sub- 
jects: exoplanet detections; solar eruptions; variable and 
pulsating stars; hi gh-energy accreting systems; ultr acom 



pact binaries (e.g. iDesidera et al.ll201ll; iFarrell et al. 



Foullon, Vcrwichtc & Nakariakov 2009; Hakala ct al 



Nataf Stanek fc Bakoj|2010l; iNess et alj|201ll;l0miva et all 



201 



200 



201 ll; iQian Solomon fc Mlynczakl l201oi ISartv et all 1200a 
Uthas et al.ll2012l : IXu et allhoilUWen et al.ll2006f h 



Figure [3] presents Lomb-Scargle plots of four stages of 
the light-curve: during the early decline (4 x 10 5 < t < 
4.7 x 10 6 s), late decline (4.7 x 10 6 < t < 9 x 10 6 s), the 
plateau (9 x 10 6 < t < 1.25 x 10 7 s), and the recent post- 
plateau stage (t > 1.25 x 10 7 s). We used a freely available 
IDLcodeS i mplementing the formula tion of lPress fc Rvbickil 
(| 19891 ) and[H orne fc Baliunasl (| 19861 ). Power and detection 
thresholds assume a grid of 1000 evenly spaced frequencies. 
For the sake of smoother curves, the curve is drawn with 10 4 
intermediate frequencies (not involved in the calculation of 
detection thresholds). In our plots, peak detection thresh- 
olds are marked in blue and red. Higher peaks are statis- 
tically significant at the 1% level ("false alert probability" 
FAP= 0.01). The blue (dashed) thre shold is zq in equation 
(18) of section III(c) of lScargld (jl982l) . The red (dot-dashed) 
line is the threshold obtained from 10 4 "white noise simula- 
tions" conservatively taking the light-curve's total variance. 
The data are not de-trended, so the variance is an overes- 
timate. The uncorrected decay timescale (to ~ 2 x 10 6 s) 
during the early and later decline (t < 9 x 10 6 s) is prob- 
ably responsible for artefacts at the long-period end of the 
periodograms (top panels of Figure [3|. 

At early times (t < 4.7 x 10 6 s) there are peaks of var- 
ied significance at periods of r as 0.23, 0.45, 0.9 and 1.4 Ms. 
The first two in this series correspond to the two features 
tentatively identified in the Fourier spectrum (§ 13.2. 1 [) . A 
characteristic spacing « 0.45Ms between dips is also in- 
ferred from a visual inspection of the early part of the 
XRT lightcurve (Figure]?]). During later weeks of the de- 
cline (4.7 x 10 6 < t < 9 x 10 6 s) there is a stronger 
peak at period r ~ 0.9Ms. Later, during the plateau stage 
(9 x 10 6 < t < 12.5 x 10 6 s), the power of the peaks is rela- 
tively low, consistent with the visible steadiness of the light- 
curve. There may be a feature at r ~ 1.1 x 10 6 s, which is 
near the 1% significance level, and perhaps another, weaker, 
feature at ~ 0.4 x 10 6 s. In the very latest times after the 
plateau, when dipping resumed (t > 12.5 x 10 6 s), there is 
a strong peak at r ps 1.4Ms. This is consistent with a visual 
inspection of the last section of the lightcurve, as it is the 
characteristic spacing between major quasi-sinusoidal dips 
(Figureg]). 

In summary, the Lomb-Scargle analysis gives signs of 
periodicity (which are stronger in some intervals than oth- 
ers), and these indications are clearer than from the Fourier 
analysis. We note that the Lomb-Scargle periodogram is 
by construction most sensitive to oscillations similar to a 
sine wave. Repetitive patterns that are highly non-sinusoidal 
may have their true significance understated. The dips seen 
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Figure 4. Top panel: zoomed-in view of an early section of the 
light-curve, showing a structure of recurrent dips, many of them 
consistent with a quasi periodic timescale of rs 4.5 X 10 5 s (solid 
green lines). This kind of dip structure produces the characteris- 
tic frequency signal detected in the Lomb-Scargle periodograms 
and structure function analyses. However, some expected dips 
are skipped or premature (dashed green lines). Bottom panel: 
zoomed-in view of the last section of the light-curve observed to- 
date; the dips occur on a longer timescale. 



in the early phases of the XRT lightcurve (Figure |4| are 
sharp and brief compared to the non-dip conditions: this is 
clearly far from a sinusoidal pattern. This motivates further 
study with another, broader technique. 

3.2.3 Structure function technique 

To obtain a more robust test of the presence or absence of 
periodicity in the X-ray luminosity variations, we carried out 
a more effective alte rnative analysis based o n Kolmogorov's 
structure functions (|Kolmogorovlll94ll . Il99lh . They are the 
statistical moments of a temporally varying signal, with the 
light curve compared to itself at an offset r. The n th -order 
structure function is 



S»(r) = <[*(t + r) -*(*)]"> 



(1) 



The second order structure function, &, is particu- 
larly useful for the purpose of time series analysis, as 
it describes the variance of a signal on timescales t. 
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Table 1. Timcscalc of regular features detected in the S 2 
structure function for data extracted in the intervals Ati = 
[0.4.8.9] X 10 6 s and At 2 = [1.0,8.9] X 10 6 s. Minima are tab- 
ulated if their depths exceed la and 3a locally, and are not en- 
closed within the catchment of another local minimum. To em- 
phasise the periodicity, we express the timescales r as multiples 
of a characteristic to = 0.445 X 10 6 s. 



In addition to their original applications in t he physics 
of turbulence |Brandenburg fc Nordlundl l201ll . and ref- 
erences therein), structure functions have a venerable 
history in the diagnosis of variability in blazars, other 



(e.g. Simonetti, Cordes & Heeschcn 1985; Cordes & Downs 


19851; Bregman et all 1988 


; Hughes, Aller & Ailed 


1992; 


Brinkmann et al. 200fJ. 


2001; Kataoka et al. 


2001; 


Ivomoto & Makishimal 200ll; Tanihata et al. 


2003; 


Wilhite et all 20081; Emmanoulopoulos. McHardv & Uttlev 


2010; Yusef-Zadeh et al. 201 ll; Voevodkin 201 ll). In 


these 



contexts they are most often used to search for break fre- 
quencies that are su pposed to characterise internal scales of 
the sy stem, though lEmmanoulopoulos. McHardv fc Uttlevl 
|201Cf ) warn users to beware of potential artefacts from 
sparsely sampled data. (Our data are densely sampled, 
even during the dips, and we are seeking signs of repetition 
rather than break scales.) Structure functions are also used 
to diagnose simulations of turbulence in interstellar media 



llBoIdyrev. Nordlund fc Padoanl |2002|; iKritsuk et al.l |2007t 
120091 ; ISchmidt. Federrath fc Klessenl l200&t). and syntheti 



obser vables of simulated jets and AGN ( Saxton et al.ll2002l . 
bOld ). In the Appendix, we show the procedures by which 
we compute the structure functions of various orders from 
data sets with unevenly sampled data. We also show how 
the corresponding uncertainties are determined. In the rest 
of this section, we report our main findings. 
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Figure 5. Structure functions of the Swift/XRT count rate at or- 
der n = 2, 4, 6, 8 (red, yellow, green and blue curves, respectively) 
for the light-curve interval 4 X 10 5 s < t < 4.7x 10 6 s. The thickness 
of each curve indicates la uncertainties. Top panel: result based 
on the direct XRT count rate, without correcting for the long- 
term decline. We interpret the rising trend towards the longest 
timescales (r > 10 6 s) as an artefact arising from the long-term 
dimming in the X-ray luminosity. Middle panel; result when the 
count rate is de-trended. Locally there is an e-folding timescale 
of 2 X 10 6 s. Bottom panel: same as the middle panel, but plotted 
in logarithmic scale, to highlight the essentially featureless "red 
noise" at variability timescales r < 10 s s. For r < 10 s s, but 
longer than the orbital period of Swift 5800 s), the structure 
functions show a featureless power-law form S n oc r n / 2 , indicat- 
ing stochastic red noise. 
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Figure 6. Top panel: structure functions of the de-trended X- 
ray count rate at order n = 2,4, 6, 8 (red, yellow, green and blue 
curves, respectively) during the time interval 4.7 X 10 6 s < t < 
9.0 X 10 6 s. A period at r « 9 X 10 5 s (and its multiples) is 
the dominant feature. The shorter periods that dominated ear- 
lier stages of the light-curve (Figure [5} arc only slightly visible 
now, and only in higher-order structure functions. Bottom panel: 
structure functions in logarithmic scale, showing the white and 
red noise at small values of r. 



As the X-ray light-curve is not a simple exponential or 
power-law decline, but seems to go through phases of differ- 
ent slope, we calculated separate structure functions during 
different time intervals, to determine whether the charac- 
teristic variability timescale changed in different epochs. We 
calculated two sets of functions: one with the count rates di- 
rectly detected by the XRT; the other with a de-trended 
count rate, to remove the long-term evolution (timescale 
of weeks/months) and highlight the short-term variability 
(timescale of a few days). At early epochs (in particular, 
at times 2.0 X 10 5 s < t < 9.0 x 10 6 s; Figure 1), multi- 
plying the count rate by an exponential function exp(t/to), 
with to ~ 2.0 x 10 6 s, is a simple way to correct for the long- 
term evolution. For the whole light-curve (including the later 
plateau stage), a simple exponential or power-law function is 
not applicable; we used empirically fitted rational functions. 
Note that the results that we obtained are independent of 
the precise choice of normalising function for the count rate, 
and including or omitting data from the first few days also 
does not significantly alter the outcome. Our choice of bin- 
ning (for example, using the same binning of the snapshot 
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Figure 7. As in Figure^ for 9.0 X 10 6 s < t < 12.5 X 10 6 s. 

observations, or a fixed time binning of 100 or 200 s) does 
not have significant effects, either. 

What we are looking for in a structure function 5„(r) 
are local depressions around specific r values: they indicate 
that the light-curve is repetitive on that period. We obtain 
best estimates of the periods by numerically locating the 
minima in 52 (the smoothest structure function) and then 
discarding all those that are shallower than the local un- 
certainty IAS2 (lcr result) or 3AS2 (3a result). We discard 
the lesser minima that are within the la (3a) catchments 
of wider and deeper minima. Table [1] shows results obtained 
from data spanning the decaying phase of the light-curve. 

3.2.4 Variability during the decay phase 

The most significant result of our analysis is that for 
t < 9 x 10 6 s, all structure functions exhibit strong min- 
ima at integer and half-integer multiples of a basic period 
r ~ 4.5 x 10 5 s (see Figure [5] for the very early epoch, 
and Figure [6] for a later epoch) . Most of these features are 
many times deeper than the sizes of the statistical uncer- 
tainties. Corresponding features appear in structure func- 
tions of different order. There is also a significant differ- 
ence between structure functions calculated for early epochs 
(t < 4.7 x 10 6 s) and later in the decline phase (4.7 x 10 6 s 
< t < 9.0 x 10 6 s). At early epochs, the most prominent 
feature is in fact at r w 2.3 x 10 5 s, that is half of what we 
identified as some kind of fundamental timescale. Figure [5] 
shows 17 consecutive features in the structure functions, at 
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Figure 8. Structure functions of the detrended X-ray light-curve, 
for the recent epoch following the end of the plateau stage. The 
top panel is integrated over the whole available sequence, (t > 
1.25 X 10 7 s). The bottom panel uses only a shorter interval (1.25 X 
10 7 < t < 18.5 X 10 7 s) before the observations became too sparse. 

integer multiples of 2.3 X 10 5 s. At later times, the shorter 
periods (r « 2.3 x 10 5 s and r « 4.5 x 10 5 s) fade and 
remain only marginally detectable in the higher order struc- 
ture functions. Instead, the deepest features in the structure 
functions are at r « (0.9, 1.8, 2.7, 3.6) x 10 6 s, higher multi- 
ples of the fundamental timescale (Figure [6}. 

The characteristic timescales found from structure func- 
tion analysis (in particular, the features at 2.3 x 10 s, 
4.5 x 10 5 s and 9 x 10 5 s) confirm and strengthen the re- 
sults obtained from the Lomb-Scargle periodogram (§ l3.2.2p 
and are consistent with a visual inspection of the lightcurve 
(Figure [4j. The two shortest frequencies correspond also to 
the two features we identified as marginally significant in the 
power density spectrum (§ 13.2.11) . None of these timescales 
appear to be exact periodicities: sometimes a major dip may 
commence slightly earlier or later than this time interval, or 
may not appear at all, but subsequent dips tend to return 
to the original phase (see intervals marked in the top panel 
of Figured]). 

3.2.5 Variability during the plateau phase 

After t « 9 x 10 6 s, the slope of the long-term trend flat- 
tened and the light-curve settled on a plateau (Figure [T]), 
with an X-ray count rate « 0.1 ct s _1 . At the same time, 



major dips, which had punctuated the light-curve until then, 
appeared to vanish. Structure functions calculated for the 
plateau stage do not show any significant periodicity (Fig- 
ure [7|). Depending on the specific choice of detrending func- 
tion, the fluctuations are "white noise" (flat shape) or flat- 
tish "red noise" (rising at long r) for timescales r > 5 x 10 4 s. 
As always, at shorter timescales the slopes of the structure 
functions indicates red noise down to the orbital period of 
Swift and the typical gap interval between observations. 



3.2.6 Variability during later epochs 

A dipping behaviour began again at t « 1.25 x 10 7 s. The 
range of count-rate variability is about the same as in the 
early phases, that is a factor ~ 10-20 between peaks and 
troughs, but the dips are broader. The dip duration has 
become comparable to the inter-dip duration. In the early 
epochs, the typical duration of a dip is ~ 1 day; at later 
epochs, it is ~ 10 d (cf. top and bottom panel of Figure Q. 

Figure [8] shows the structure functions of the detrended 
light-curve during the most recent epoch. There are minor 
wiggles at O.lMs timescales. These are probably spurious, 
as they are comparable to the gaps between observations. 
When the whole time interval t > 1.25 x 10 7 s is included in 
the structure function calculation (top panel of Figure [S]), 
it is not immediately obvious whether there are dominant 
timescales. This is partly because the Swift monitoring ob- 
servations have become less frequent, and partly because 
the source has not varied much after t w 1.85 x 10 7 s. Tak- 
ing only the sub-interval 1.25 x 10 7 < t < 1.85 x 10 7 s, 
when the source is more variable, reveals clearer features at 
r w 0.7,1.4,2.9,4.3,5.0 and 5.7 Ms (bottom panel of Fig- 
ure [SJ — all simple multiples of the same basic frequency, 
although it is a different base frequency than the one as- 
sociated with the narrow dips in the initial rapidly declin- 
ing stages. The variability timescale r ~ 1.4 x 10 6 s was 
also identified at high significance in the Lomb-Scargle peri- 
odograms (fourth panel of Figure[3]). In fact, it appears more 
significant in the Lomb-Scargle periodograms than in the 
structure functions. This is probably because of the nearly 
sinusoidal waveform of the recent oscillations. 



3.3 Time-dependent spectral properties 

3.3.1 Hardness ratio evolution 

It was noted in previous work (|Burrows et al. 

I 1201 ll : 

iLevan et alll2011al ) that the X-ray emission of Sw J1644+57 
is harder (flatter spectrum) when the source is brighter. 
However, this one-dimensional correlation cannot fully char- 
acterise the spectral behaviour. The count rate drops during 
the short-term dips, but also as a result of the long-term evo- 
lution: we need to treat the two effects separately, as they 
probably have different physical origins. To do so, we di- 
vided the light-curve into several time intervals, similarly 
to our approach with the Lomb-Scargle periodograms and 
structure function analysis; the intervals are also defined in 
approximately the same way for convenience, although the 
precise start and end of each phase is somewhat arbitrary. 
We plotted the hardness ratio (defined as the observed count 
rate in the 1.5-10 keV band over that in the 0.3-1.5 keV 
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band) as a function of total count rate for each segment, in 
different colours, in Figure [5] 

At early times, the general trend is a moderate hard- 
ening as the baseline flux declines, alternating with signifi- 
cant softening during each dipping episode. This is partic- 
ularly evident in the distribution of datapoints coded with 
red, green and magenta colours in the top two panels of 
Figure [9] which cover the time interval 2 x 10 s < t < 
4.7 x 10 6 s. Further independent confirmation of the soften- 
ing behaviour during the dips comes from XMM-Newton ob- 
servations taken on 2011 April 16, April 30, May 16 and May 
30. The power-law photon indices measured during the first, 
second and fourth observation were V m 1.88, 1.66, 1.61 re- 
spectively (|Miller fc Strohmaver|[2l)l"lh . consistent with the 
moderate long-term hardening as the baseline flux declined; 
but the third observation occurred during a dip, and the 
photon index was F « 1.97. 

At later times, as the light-curve reaches a plateau and 
dips become less frequent, the hardness ratio also seems to 
settle around a constant value corresponding to T ~ 1.4, 
with some scatter mostly due to the small number of counts 
in each snapshot observation. When the dipping behaviour 
resumes, at time t « 1.25 x 10 7 s, the hardness ratio no 
longer changes during the dips (bottom panel of Figure [9]), 
in contrast with the behaviour at early times. 

In summary, we found that at early epochs, the dips 
are shorter, sharper and significantly softer than the base- 
line emission at the same epoch. At late times, instead, the 
dips become broader and longer, and the hardness ratio is 
independent of flux. We sketch our interpretation of the flux- 
hardness evolution in Figure [TD] 



3.3.2 Spectral evolution 

Our investigation of hardness evolution makes it clear that 
we need to distinguish between the short-term and long- 
term spectral evolution. To do so, we use the same method 
applied to our structure function analysis: we normalised the 
snapshot light-curve to a smooth empirical function that 
fits the long-term evolution. We then extracted intensity- 
selected spectra based on the normalised count rates rather 
than the observed ones; for example, integrated spectra from 
all time intervals in which the normalised flux was > 0.9, 
or between 2/3 and 0.9, etc. We also fixed the flux level 
and extracted time-selected spectra, that is with the same 
normalised flux band but at different epochs. 

In Table [2] we summarise the result of spectral fitting 
for a sample of intensity-selected spectra from the decline 
phase (t < 9 x 10 s). All spectra are well fitted by ab- 
sorbed power-laws, with a column density ~ 1.5 x 10 22 cm -2 
in the local frame at redshift 0.354. Small variations in 
the fitted column density from band to band may not be 
statistically significant, and do not appear to follow any 
trend. The high- and medium-intensity bands have a simi- 
lar photon index ~ 1.7. Low- intensity bands have a steeper 
slope, as already suggested by our hardness ratio analysis: 
as steep as T as 2.9 for spectra taken across the bottom of 
the dips (Figure lll[) . This result confirms our interpreta- 
tion of the hardness- intensity plots and is consi stent with 
the spectral results of lMiller fc Strohmaverl |201 1[ ). It shows 
even more clearly that the early dips are characterised by 
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Figure 9. Top panel: hardness-intensity diagram for the snap- 
shot datapoints in PC mode. Red datapoints are those for t < 
1.25 X 10 6 s; green datapoints for 3.41 X 10 6 < t < 4.7 X 10 6 s; 
blue datapoints for 4.7 X 10 6 < t < 9x 10 6 s; orange datapoints for 
9 X 10 6 < t < 1.25 X 10 7 s. Error bars have been omitted for clar- 
ity. The horizontal dashed lines correspond to the hardness ratio 
expected in PC mode for a power-law spectrum with intrinsic col- 
umn density N K = 1.5 X 10 22 cm -2 (at redshift 0.354), a Galactic 
absorption = 2 X 10 20 cm -2 , and several different values of 
the photon index T. Middle panel: hardness-intensity diagram for 
the datapoints in WT model, at time 1.25 X 10 6 < t < 3.41 X 10 6 s, 
with intrinsic ./Vjj = 1.5 X 10 22 cm -2 . Bottom panel: hardness- 
intensity diagram in PC mode for the last part of the light-curve 
(1.25 X 10~ < t < 2.1 X 10 7 s), binned to observation datapoints; 
an intrinsic TVjj = 1.9 x 10 22 cm -2 was used for this late epoch, 
as suggested by our spectral fitting (sec § 13. 3.211 . 
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Parameter 


Band 1 Value 


Band 2 Value 


Band 3 Value 


Band 4 Value 


Band 5 Value 


Model: tbabs*ztbabs*power- law 


^Vh, 2=0. 354 

r 


[0.02] 
i k.7+0.10 

1 -°'-o.io 

1 67+ ' 06 
16.61! .2 


[0.02] 
1 41+O.09 

, 7 9 +0.06 
1 - ' -0.06 

1 n n + 1-0 
12 -°-0.9 


[0.02] 

1 qc+0.06 
1 -°' J -0.06 

1 KQ+004 
i - DM -0.04 

O+0.4 
8 - 2 -0.4 


[0.02] 

1 09+0.14 
J -° -0.12 

, 09+O.O9 
1 -° z -0.09 

„ c+0.4 
o-C-0.4 


[0.02] 
1 46+ -44 

9 07+O.47 
^■»'_0.40 

1 c + 1.1 
i -8-0.6 


/0.3-10 
-£-0.3-10 


7 7+O.I 
'•'-0.2 

46.8+J' 3 


■5 o+O.l 


3-8t ;i 

99 c+0.6 


1 o+O.l 
J - a -0.1 

a q+0-6 
8-9-0.5 


20+ 02 

r O+4.0 
b.8_2. 8 


xl 


1.06 

(328.3/311) 


0.96 
(271.1/282) 


0.96 
(461.1/480) 


0.98 
(155.2/158) 


1.03 

(26.8/26) 



Table 2. Best-fitting spectral parameters for de-trended intensity-selected Swift/XRT-PC spectra. Band 1 is for time intervals when the 
count rate was > 9/10 of the spline-fitted baseline rate (Figure[T}; band 2 for count rates between 2/3 and 9/10; band 3 for count rates 
between 1/3 and 2/3; band 4 for count rates between 1/10 and 1/3; band 5 for count rates < 1/10 of the baseline. All bands include only 
intervals at 2 X 10 5 s<i<9x 10 6 s. Units: Ah, Gal an d A r H , z =o.354 are m units of 10 22 cm" 2 ; N po in 10~ 3 photons keV" 1 cm~ 2 s _1 
at 1 keV; fo.3— 10 in 10 — 11 erg cm~ 2 s ; Lo.3-10 in 10 45 erg s -1 . Errors indicate the 90% confidence interval for each parameter of 
interest. 



Parameter 


Band A Value 


Band B Value 


Band C Value 


Band D Value 


Model: tbabs*ztbabs*power- law 


A H ,Gal 
-^H,z=0.354 

r 


[0.02] 

1 50+ ' 07 
i - ou -0.09 

-, 77+O.O5 
L -< '-0.05 

30.2±!;| 


[0.02] 

1 54+0.19 
l- o *-0.20 

1 41+0.11 

1 6+°- 3 
l-°-0.2 


[0.02] 

1 79 +0.06 
l-' 3 -0.09 

l- 4u -0.04 

1 o+O.l 
1- J -0.1 


[0.02] 

1 Q9+0.39 
!- MZ -0.34 

l- ou -0.17 

1 0+°' 3 
l- u -0.2 


/0.3-10 

£-0.3-10 


12 9+0.1 

lz - z -0.2 

so-sl 2 ,^ 


1 1+°1 

i-i-o.i 

e 7+0-3 
°-'-0.3 


n 01 +0-05 
u - y l-0.05 

4 9+0-1 


o4+0- 04 

9+O.2 
°- -0.2 


xl 


1.06 
(432.3/408) 


0.90 
(95.6/106) 


0.92 
(487.5/528) 


1.06 
(78.4/74) 



Table 3. Best-fitting spectral parameters for time-selected Swift/XRT-PC spectra. All spectra include only intervals in which the 
count rate was > 2/3 of the normalized continuum rate. Band A is for time intervals 2 X 10 5 s < f < 4.7 X 10 6 s; band B for 
4.7 X 10 6 s < t < 9 X 10 6 s; band C for 9 X 10 6 s < t < 1.25 X 10 7 s; band D for 1.25 X 10 7 s < t < 2.1 X 10 7 s; Units arc as in TableQ] 
Errors indicate the 90% confidence interval for each parameter of interest. 



a steeper spectral slope, not by higher absorption: they are 
not eclipses by orbiting or rotating clouds. 

We then considered high- and medium-intensity spec- 
tra (which, as we have seen from Table have the same 
average spectral slope), taken from all time intervals with 
normalised count rates > 2/3, and divided them into sev- 
eral time bands corresponding to the intervals used for tim- 
ing analysis. We found (Table [3| that all spectra are well 
fitted by a simple absorbed power-law; adding blackbody 
or optically thin thermal components does not improve the 
fit, nor is there any evidence of breaks in the power-law at 
high energies. The slope in the early epoch (t < 4.7 x 10 6 s) 



was r « 1.8, flattening to F « 1.4-1.5 at later epochs (Fig- 
ure 1 1 2 p . This long-term spectral trend is in perfect agree- 
ment with what was inferred from the hardness-intensity 
study (§ 13. 3. H and Figure [10)) . The intrinsic column density 
appears to increase at later epochs (Table [3}, although this 
result is only marginally significant. 



At late epochs, the unabsorbed isotropic 0.3-10 keV lu- 
minosity becomes as low as ~ 4 x 10 44 erg s _1 during the 
dips (TableQ. This is now approaching the pre-burst up- 
per limit to the X-r ay luminosity, Lx < 1.7 x 10 44 erg s _1 
|Bloom et alj|2011bj |. If Sw J1644+57 continues its decline, 
it is plausible that in a few months' time we will start see- 
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Figure 10. Schematic interpretation of the hardness-intensity 
plots. The fast decline during the first 3 months is associated to 
a moderate hardening (from T as 1.7—1.8 to T Ri 1.4) followed 
by a plateau. At early stages, the spectrum becomes much softer 
during each dip; at later times, the hardness ratio does not change 
in the dips. 



Parameter 


Value 1 


Value 2 


Model 1: tbabs*ztbabs*power-law 
Model 2: tbabs*ztbabs*raymond-smith 




[0.02] 


[0.02] 


^H,z=0.354 
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N po 


, cro+0.73 

i.oo_ Q70 
14+ 007 


1 7+ 1 - 2 
i -'-0.9 


kT c 




> 3.9 


N IS 




7 o+5.0 
'• 8 -1.9 


/o.3-10 
£-0.3-10 


n n 7 o+0.013 
U.UM_q 012 

45+ 30 


064+ 011 
U.UD^_ 032 

39+ 08 
U - OM -0.04 


Cstat 


7.35/12 


6.16/12 



Table 4. Best-fitting spectral parameters for the Swift/XRT-PC 
spectrum during the late dips (1.25 X 10 7 s < t < 2.1 X 10 7 s), 
fitted with the standard power-law model and with an optically- 
thin thermal plasma model. Units of JVh, flux and luminos- 
ity arc as in Table [2] The normalisation N IS = 10 — 11 (1 + 
z) 2 /(4-7T D 2 ) J{n a n\dV) where D is the luminosity distance. Er- 
rors indicate the 90% confidence interval for each parameter of 
interest. 



ing the unbeamed component of the emission (e.g., from an 
accretion disk or its corona), or we will infer stronger limits 
on the unbeamed flux. We tried fitting the late-epoch dip 
spectrum with a thermal-plasma model and found that it is 
at least as good as a power-law model (Table [4]); there are 
not enough counts to discriminate between them. 



4 DISCUSSION 

The general consensus is that Sw J1644+57 is a previ- 
ously quiescent nuclear BH that underwent an outburst 



caused b y accretion of a t idally disrupted star (|Bloom et al 
2011allbl: iBurrows et all l201ll : ICannizzo. Troia fc Lodato 



201ll ; Tsocratesf l201ll ). Its high apparent luminosity requires 



strong beaming. The emission in the keV band can be ex- 
plained as synchrotron and/or inverse Compton radiation 
from energetic electrons streaming along the relativistic jet 
pointing towards us. The relative contribution of the two 
processes is still debated. It is also still unclear how the 
jets are launched at the onset of the outburst, when a 
fully formed accretion disc cannot have had the time to 
develop, after the disruption event (cf. the violent mass 
transfer and jet formation in the X-ray binary Cir X-l; 
I Johnston. Fender fc Wull 19991 ). 

From the Swift/XBT data we have extracted further 
information on the time-dependent behaviour and spectral 
evolution of the source. We have decomposed the ampli- 
tude variations into a long-term evolution (an initial decline 
followed by a plateau) and recurrent dips. If these dips had 
been wholly stochastic processes, we would expect that fluc- 
tuations should be of the same order of magnitude as the 
depths of the dips. Otherwise it would require many minia- 
ture flares (a day's worth of them) to coherently conspire to 
become small simultaneously, or to abstain simultaneously. 

However we have shown that the dips are not merely 
random fluctuations (as is often the case in the X-ray light- 
curve of accreting sources), but have a certain regularity, 
apparently governed by an underlying periodic driver. We 
have used Lomb-Scargle periodograms and structure func- 
tion analysis to search for the main characteristic timescales 
in the frequency range ~ l(T 7 -l(n 4 Hz. We expect vari- 
ability on those timescales based on physical arguments. 
For example, there may be a residual inhomogeneous dis- 
tribution of debris at or near the tidal disruption radius 
Rtd w WGM./c 2 x [M./(10 7 M Q )] _2/3 (for a solar-type 
star). The Keplerian times cale £k — 2t t/Q,k at the tidal dis- 
ruption radius is ~ 10 4 s (|Rees|[l98Sl ). independent of BH 
mass in the Newtonian approximation. The ring of debris 
is likely to be still oriented in the orbital plane of the dis- 
rupted star: this may generate a Lense-Thirring precession 
with a timescale i L T = \t K /( 2a)} x \(R TD c 2 ) / '(GM.)f /2 ~ 
30-300t K /(2a) ~ 10 5 -10 6 s (lArmitage fc Nataraianl 1 19991 ; 
iMerritt et aHl20ld : IStone fc Loebll2012r ). where a is the BH 
spin parameter. We found significant features in the Lomb- 
Scargle periodograms and structure functions, especially at 
t m 2.3 Ms, 4.5 Ms, 9 Ms during the early epochs, and 
t m 0.7 Ms and 1.4 Ms during the late epochs. It appears 
that the characteristic dipping timescale (and its associated 
harmonics) has shifted between the early and late section of 
the lightcurve. 

Whatever the mechanism for the dips, it is something 
that also changes the spectral properties of the emission, as 
the spectrum is steeper during the early dips (F w 2.8-3) 
than immediately before and after (r ~ 1.4—1.7). At the 
same time, the long-term trend outside the dips shows spec- 
tral hardening at the beginning, followed by a plateau at 
T « 1.4. 

Here, we attempt to make sense of our findings in the 
framework of a synchrotron-emitting relativistic jet scenario. 
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Figure 11. Best-fitting power-law models and x 2 residuals 
for three de-trended intensity-selected spectra, showing spec- 
tral steepening in the dips. Blue datapoints and residuals: 
background-subtracted spectrum extracted from all the time in- 
tervals when the count rate was > 9/10 times the de-trended 
baseline (Band 1 in Table [2}. Green datapoints and residuals: 
spectrum from the time intervals when the count rate was be- 
tween 1/10 and 3/10 of the de-trended baseline (Band 4 in Ta- 
ble [2). Red datapoints and residuals: spectrum from the time 
intervals when the count rate was ^ 0.1 times the de-trended 
baseline (Band 5 in Table [2}. For all three spectra, we used only 
data up to t = 9 X 10 6 s; we used only data in PC mode, to reduce 
background contamination. See Table [2] for the fit parameters. 




Energy (keV) 

Figure 12. Best-fitting power-law models and x 2 residuals for 
three time-selected spectra, showing moderate long-term harden- 
ing. Red datapoints and residuals: background-subtracted spec- 
trum extracted from the time interval 2 X 10 5 s < t < 4.7 X 10 6 s 
and only including intervals with a count rate > 2/3 times the 
de-trended baseline. Green datapoints and residuals: background- 
subtracted spectrum at 4.7 X 10 6 s < t < 9 X 10 6 s and count 
rate > 2/3. Blue datapoints and residuals: background-subtracted 
spectrum at 9 X 10 6 s < t < 1.25 X 10 7 s and count rate > 2/3. 
All data are in PC mode. Sec Table [3] for the fit parameters. 



We propose that: (i) the jets have precessional and nuta- 
tional motion; (ii) the jets have a collimated core surrounded 
by an envelope of less energetic, less collimated electrons 
straying out of the core; (iii) the synchrotron-emitting jets 
become more compact as the source declines, after the initial 
outburst. 

In this scenario, the baseline flux in the XRT band is 
direct synchrotron emission from a jet pointing toward us. 
In the first few days after the outburst, when the accre- 
tion rate is highest, the jet expands and propagates for- 
ward into the interstellar medium. The X-ray emission re- 
gion is most extended at that early epoch, and it is opti- 
cally thin to synchrotron radiation. Energetic electrons are 
freshly accelerated in the jets, with an energy distribution 
N(E) oc E~ p , where the energy spectral index 2.2 < p < 2.5 
(see e.g jBallard fc Heavens]|l992l : lKirk fc Duffy||l999f ). The 
synchrotron emission spectrum from this population of elec- 
trons has a specific intensity I v oc v~ a with a ~ 0.7, cor- 
responding to a photon index T ~ 1.7. When the most vio- 
lent flarings subside and accretion becomes more steady, the 
overall accretion luminosity of the system decreases and the 
jet power is reduced accordingly. In this phase, synchrotron 
emission no longer comes from the expanding jet ejecta, but 
it is rather confined to the vicinities of shocks at the base of 
the jet. As the X-ray emission region becomes more com- 
pact and denser, the synchrotron photons start to suffer 
self-absorption, leading to the observed spectral flattening: 
approaching a similar situation to the flat-spectrum core 
emission in radio galaxies with compact jets. 

Moreover, the transition from extended to compact jet 
emission may also provide a simple explanation for the lev- 
elling off of the X-ray luminosity at late epochs. For an opti- 
cally thin source, the luminosity depends on the total volume 
of the synchrotron-emitting region in the jets; a constant 
luminosity would require a high degree of fine-tuning. By 
contrast, in an opaque source, the luminosity depends on 
the effective area of the emission region visible to us. The 
growth of the radio emitting region is perhaps explained by 
a more diffuse outer region, such as the larger jet-driven tur- 
bulent cocoon and the expanding bow shock surrounding the 
whol e system |l3ower. Bloom fc Cenkdl201ll : IZauderer et alj 
l201ll : iBurrows et alj|201ll ). Modelling the radio flux evolu- 
tion requires knowledge of the jet parameters, which are un- 
fortunately not well constrained by the observational data 
obtained so far. We therefore leave the jet modelling, which 
is beyond the scope of this paper, to a future study. 

Regarding the dips, we showed that they are not due to 
partial occultation of the jets by some line-of-sight material, 
as we would expect higher absorption rather than spectral 
steepening in this situation, contrary to what is observed. 
Instead we propose that dips occur when the synchrotron- 
emitting jet is not perfectly aligned along our line-of-sight, 
and hence the most highly beamed component of the syn- 
chrotron emission does not reach us. The observed X-ray 
flux is then dominated by the lower level of emission from 
the slower and/or less collimated electrons in a sheath or 
cocoon around the jet. This emission could be synchrotron 
from an aged (steeper-spectrum) population of electrons, if 
there are significant magnetic fields in the jet envelope, or 
synchrotron self-Compton. We must also stress that even at 
the bottom of the dips, the apparent luminosity is higher or 
similar (for the late dips) to the Eddington luminosity of a 



© 2011 RAS, MNRAS OOCHHTrTl 



X-ray variability of Swift J1644+57 13 



10 7 Mq BH. This suggests that the emission is moderately 
beamed even during the dips. 

In addition to the expected long-term decline and 
stochastic flux variability, it is intriguing that many of the 
dips occur on regular patterns, likely to be governed by a 
periodic driver; it is even more intriguing that sometimes a 
dip is skipped but the following one happens again approxi- 
mately at the expected phase (Figure |¥J. Careful inspection 
of the centroid locations of the dips has revealed that they 
are not exactly phase-matched but can be slightly advanced 
or delayed. To explain these peculiar dip patterns, we pro- 
pose that the jet undergoes nutation as well as precession. 
The combination of the two effects drives the cone of the jet 
out of our line of sight, from time to time. The presence of 
nutation is essential, as it naturally explains why the dips 
occur in a regular but not exactly periodic pattern, and why 
there are slight advances or delays of the dip centroids. 

We note also that the appearance of the dips changes: in 
the early epoch we see a slow fainting and rapid rise; but in 
later dips the dimming and rising are symmetric (Figure [4}. 
The fractional time that the object is in a dip or out of 
a dip also changes: brief dips at early times; long dips in 
late epochs. These observations gain significance in the jet 
precession/nutation scenario, especially if the emission core 
became more compact. When the early dips were sharp and 
brief, the precession subtended a tight angle, compared to 
the jet's beaming angle. The larger width of the dips at 
late times, and the larger fraction of time spent in a dip, 
probably indicates that the core of the jet is narrower (or 
swinging more widely) while the layer of slower electrons 
around the core grows in size. The asymmetric dipping at 
early times may mean that the region around the core is not 
moving in phase with the jet core; lagging while the core 
moves out of the line of sight, while being compacted when 
the core moves back. 

The remaining question is what causes the jet precession 
and nutation. Jet precession is expected as the BH spin axis 
is unlikely to be perfectly aligned with the initial angular 
momentum of the disrupted star. Thus, the normal direc- 
tion to the plane of the transient accretion ring/disc formed 
by the stellar debris is tilted with respect to the BH spin, 
and this causes the jet to precess. In addition, jet nutation 
can be triggered if the precessing jet/disc system has been 
perturbed, either by an internal or external driver. In the in- 
ternally driven case, nutation may be induced by an uneven 
mass distribution in the transient debris ring/disc. In the 
externally driven case, the jet and the accretion ring/disc 
may be perturbed by the orbiting remnant of the partially 
disrupted star — for example, a white dwarf on an ellipti- 
cal orbit would make several passes inside the tidal radiu s 
before being completely disrupted (|Krolik fc Piranl l201ll ). 
Alternatively, they may be perturbed by another gravitat- 
ing object, such as (unlikely but not impossible) a satel- 
lite massive black hole (as predicted in nuclear swarms of 
compact objects: Freitag. Amaro-Seoane fe Kalogeral 120061 ; 



iPortegies Zwart et al.ll2006l ): or perhaps a more supermas 
sive companion (about which Sw J1644+57 is the satellite). 

An important clue to the evolution of the system 
comes from the change in the characteristic variability 
timescales at early and late times. If the nutation is ex- 
ternally driven, the change of period might conceivably be 
explained if the torque changed (perhaps during an ec- 



centric perimelasma passage). Scenarios dominated by disc 
torques suggest other possibilities. A possible explanation is 
that the warped disk has not yet reached a steady state 
or Bardeen-Petterson r egime (|Bardeen fc Pettersonl Il975l : 
Kumar fc Pringldll985h . In other words, the warp radius 
(transition between the inner part of the disc aligned with 
the BH spin and the outer part aligned with the orbital plane 
of the mass donor) may still be propagating outwards, and 
as a result the Lense-Thirring precession timescale at the 
warp radius is incr e asing. Alternatively, it was suggested 
(iFragile et al.l 120071 ; iDexter fc Fragile! 1201 ll : IStone fc Lo"ebl 
l2012f ) that a geometrically thick disc may never settle in 
a Bardeen-Petterson regime, and may instead precess as a 
solid body (a tilted, rather than warped disc). In this case, 
the precession frequency depends on the dimensionless ra- 
dial surface density profile of the disc, which is also likely to 
evolve in time. 

In fact, a problem of the disc-driven precessing- jet sce- 
nario is that we expect the jet axis to have already precessed 
out of alignment with our line of sight (and not just for short 
dips) after few weeks, unless the initial alignment of disc an d 
BH spin axis was already very close (|Stone fc Loebl [20121 ). 
This seems an unlikely coincidence if the disc was formed 
after a tidal disruption event. Thus, we suggest a note of 
caution regarding the tidal disruption scenario as the origin 
of the flare. It is useful to remember that the 2011 March 
28 flare marks the moment when the jet turned on, not the 
time of the putative tidal disruption event, which may have 
happened unobserved several weeks earlier, depending on 
the time required to spread the debris onto a disc, build up 
a magnetic field and launch a jet. The main argument for a 
tidal disruption as the origin of the whole process is that the 
host galaxy was not known as an AGN before the flare. How - 
ever, pre-event luminosity upper limits (|Bloom et al.ll2011bT ) 
are only implying that the bolometric luminosity was below 
Eddington for a 1O 7 M BH, or below 10 Eddington for a 
1O 6 M BH. Thus, we cannot exclude that an accretion disc 
was already present and approximately aligned, and the jet 
ejection was due to a state transition rather than a tidal dis- 
ruption. State transitions with relativistic blob ejections and 
radio/X-ray flares are common in Galactic BH transients, 
especially when they switch from the hard to the soft (ther- 
mal) state (|Fender. Belloni fc Galioll2004l ; [Fender fc Bellonil 
l2004h . Scaled with their respective mass ranges, the dura- 
tion of the flare in Sw J1644+57 so far corresponds to only 
~ 10-100 s in a typical Galactic BH. Further monitoring of 
the flux evolution over the next few months will help testing 
between the tidal disruption and state transition scenarios. 

The scenario outlined above does not imply that ev- 
ery dip is due to oscillations and precession of the jet. Our 
most general interpretation is that dips correspond to phases 
when we are seeing slower and less collimated electrons along 
our line of sight in the jet. Even when the nozzle is steady, 
jets are subject to internal instabilities of magneto-hydro- 
dynamical nature. Internal shocks echo and self-intersect 
up and down the jet, producing standing patterns with a 
characteristic spacing that depends on the Mach number 
of the flow and the width of the jet. Jets also suffer an 
axisymmetric pinch instability, a non-asxisymmetric helical 
instabilities, a nd other higher-order ha r monics of the heli- 
cal instability |llardee fc Norman] 1 19881 ; iNorman fc Hard"eel 
1 19881 ; lHardee et al.lll992TT Observed changes in the jet col- 
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limation and luminosity will depend on the interplay of the 
external driving oscillations {e.g., jet precession) and the 
resonant frequencies of the internal instabilities. The sit- 
uation becomes even more complex if — as seems to be 
the case in Sw J1644+57 (|Berger et al.ll201lT ) — the jet 
has a broad Lorentz factor distribution rather than a sin- 
gle speed. With more arduously detailed modelling work (in 
the future), tested by long-term X-ray and radio monitor- 
ing |Berger et al.ll201lT ). one might hope to infer something 
about the internal layering of the jet, and distinguish be- 
tween internal and external driven jet variability. 



5 CONCLUSION 

We investigated the luminosity and spectral evolution of 
the peculiar X-ray source Sw J1644+57 using Swift/XKT 
data. Our structure function analysis showed that the large- 
amplitude variations in the light curve of the source are 
not simply stochastic fluctuations. In particular, the oc- 
currence of at least some of the dips appears to follow 
some regular patterns characterised by multiples and frac- 
tions of a time interval r » 4.5 x 10 5 s; visually, several 
prominent dips are spaced at ~ 0.5r ~ 2.2 x 10 s s, and 
~ 2r ~ 9 x 10 J s. After the plateau the base period switched 
to t w 1.4 x 10 6 s, but the scenario remains essentially the 
same. The X-ray spectrum outside the dips is always con- 
sistent with an absorbed power-law, but its photon index 
evolves on a weekly /monthly timescale, from F ~ 1.8 at the 
beginning of the decline phase to T ~ 1.4 at a later stage 
when the luminosity seems to level off. The X-ray spectrum 
is much softer during the dips, while there is no significant 
change in the absorption column density. 

We proposed a scenario in which the synchrotron emit- 
ting jets, launched from a massive BH accreting a disrupted 
star, undergo both precession and nutation. The baseline 
flux in the keV band is direct synchrotron emission when 
the jet cone is in our line of sight; dips occur when the jet 
cone goes briefly and partially out of our line of sight. We 
argued that the synchrotron photons come from the opti- 
cally thin jet ejecta during the initial high-luminosity, ex- 
ponential decay phase; as the accretion power subsides, the 
dominant emission region becomes an opaque feature that 
is confined to a compact jet core. We attributed the fainter, 
steep-spectrum dip emission to a population of less ener- 
getic electrons, perhaps streaming in a cocoon surrounding 
the collimated jet core. Jet nutation and precession provide 
a natural explanation for the dip patterns in the X-ray light- 
curve: in particular, their preferred occurrence at some regu- 
lar intervals, their occasional disappearance and subsequent 
reappearance, and the phase advance and delay of the dip 
centroids. In addition, internal jet instabilities can produce 
oscillations in the speed and cross section of the jet flow, and 
therefore a recurrent dipping behaviour. Long-term mon- 
itoring of the characteristic oscillation frequencies will be 
required to test whether they are due to internal or external 
drivers. 
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APPENDIX 

A.l Piecewise linear light curves 

Consider a temporally evolving signal z (e.g. X-ray flux or 
counts), for which the data occur in bins i spanning tempo- 
ral intervals [U, Ui] where tj is the starting time and Ui is the 
ending time of bin i. In characterising the temporal variabil- 
ity of the light-curve, we may offset a copy of these bins by a 
lag time r, and compare to the unoffset data. With patchy, 
piecewise data, it is necessary to compute the temporal over- 
lap between each pair of unoffset and offset bins. If the time 
intervals of the offset bins are [ij,Mj] = [tj + T,Uj + r], 
then the interval of overlap with an unoffset bin [t»,Uj] is 
[Tij, Uij] = [max(tj,tj-),min(ui, u'j)]. If Uij < Ty, then there 
is no overlap. The weight or duration of overlap between a 
pair of bins i,j is 



Uij = Wij (r) = 







otherwise. 



(2) 



Suppose that the amplitudes within each bin evolve lin- 
early in time: z = ait + bi and z = ajt + bj . The linear co- 
efficients depend on the (observed) amplitudes z and times 
at the ends of each time interval, aj = (Szj)/(uj — tj) and 
bj = Zj — a jtj. During the overlap, the difference between 
the amplitudes of offset and unoffset light curves is y = 
z'jit+rj-Ziit) = AijX + Bij, with A e= (t-Tij)/(Uij -T^), 
and ^ A ^ 1 by construction, where we abbreviate 

Aij = (v'ij - v'ij)5zj - (u^ - ii t j)8zi ; (3) 

Bij = n'ijSzj + bj + ajtj — fiijSzi — bi — ciiti . (4) 

For brevity, we define relative fractions /iy , Vij , , and v[j 
expressing how far the overlap occurs along the segments i 
and j. 



Tij t>% 

Ui — U 



t, 



Tj + r 



Uij + t — tj 



(5) 



A2. Structure functions for unevenly sampled data 

For continuous data, the order-n structure function is de- 
fined as 

S n (r) EE ([z(t + T)-z(tT) 

= ±J [z(t + T)-z(t)] n dt . (6) 

© 2011 RAS, MNRAS 000,[]jfT6] 



X-ray variability of Swift J1644+57 15 



For patchy, discretely binned data, we need to integrate con- 
tributions from each pair of potentially overlapping bins, 
piece by piece. If bin i overlaps with the offset copy of bin j 
then this pair contributes an amount 7y to the numerator 
of the structure function, 



hi = hj(r) = / [z'(t + r)-z(t)] n dt 

rV tj 



[cij t + cij t + bj — ait — bi] n dt , 



(7) 



so that S n = S n (T~) = W^ 1 Iij. The normalisation fac- 
tor W is the total duration of temporal overlaps between all 
pairs of bins. 



(8) 



For m 7^ aj (and Aij 0), the integral can be evaluated 
directly, 



(n + l)Ai 



\n + l + 1 



(9) 



An alternative, binomial expansion proves to be more nu- 
merically stable in practice, especially where |Ay| is small: 



E 



/A \',i>':, '"' 
fc J /c + 1 



(10) 



A3. Determination of the uncertainties 

If each bin k has a measurement Zk, with uncertainty Azk, 
then we obtain the uncertainty on the overall structure func- 
tion by quadrature. Now to propagate the uncertainties in 
the measurements zu to uncertainties in the structure func- 
tions S n (r), one needs partial derivatives of each piece Iij 
with respect to each observable Zk- It emerges that 



dlij / n 



k Ak—l pn — fc OAjj 

fc + 1 ij ■ dz k 



n-k . h n _ k - 1 dB iJ 
+ - — — AnBi 



fc + 1 



dz k 



(11) 



If we assume a zigzag data patching scheme {z i+ \ — zt+Szi) 
then the revelant partial derivatives of Aij and Bij are 



dA l3 
dz k 

dBjj 
dz k 



{vij - iXij)5 Ki+1 + {v t j - ixij)5 k i 
VijSkj+i + (1 - v'ij)&k,j 

where 8 is the Kronecker delta symbol. Now abbreviate 
fc 



A 



a k — 1 rjn — k 

k I k + 1 XJ 13 



(12) 
(13) 

(14) 
(15) 



and 

Ci i3 = Atj (y'ij - fi'ij) + B%j n'ij ; (16) 

C 2ij = -Aij(i/ij - n'ij) + Bij(l - fi'ij) ; (17) 

C 3 y = -Aij(vij - mj) - Bijfitj ; (18) 

C-i.j = Aj^ij- -My) +Bij([iij - 1) • (19) 

It follows that 
r)l 

^ = C Uj S k j+1 + C 2tj 5 kJ + C St] 5 k:i+1 + C itj 5 fc „ . (20) 

By contracting the Kronecker deltas, we obtain 

dS„ _ 1 ^ dlij 
~dz~^~W ^ diT 



= w E ( c ' i I ,'=-i + ^t,* + C3 h-i,i + c *k,i 



(21) 



The four matrices (Ci, C2, C3, C4) only need to be computed 
once for each choice of r. If the uncertainties in the flux 
measurements are Azfc, the total uncertainty in the structure 
function is then given by 

(^) 2 -E(^ fc ) 2 . (22) 
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